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Stability analysis of chromo-natural inflation and possible evasion of Lyth's bound 
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We perform the complete stability study of the model of chromo-natural inflation (Adshead and 
Wyman '12), where, due to its coupling to a SU(2) vector, a pseudo-scalar inflaton \ slowly rolls 
on a steep potential. As a typical example, one can consider an axion with a sub-Planckian decay 
constant /. The model was recently studied (Dimastrogiovanni, Fasiello, and Tolley T2) in the 
m g >> H limit, where m g is the mass of the fluctuations of the vector field, and H the Hubble 
rate. We show that the inflationary solution is stable for m g > 2H, while it otherwise experiences a 
strong instability due to scalar perturbations in the sub-horizon regime. The tensor perturbations 
are instead enhanced at large m g , while the vector ones remain perturbatively small. Depending 
on the parameters, this model can give a chiral gravity wave signal that can be detected in ongoing 
or forthcoming CMB experiments. This detection can occur even if, during inflation, the inflaton 
spans an interval of size Ax = O (/) which is some orders of magnitude below the Planck scale, 
evading a well known bound that holds for a free inflaton (Lyth '97). The spectral tilt of the scalar 
perturbations typically decreases with decreasing m g . Therefore the simultaneous requirements 
of stability, sufficiently small tensor-to-scalar ratio, and sufficiently flat scalar spectrum can pose 
nontrivial bounds on the parameters of the model. 



I. INTRODUCTION 



Inflation is a successful paradigm for the physics of 
the early Universe [l[ . Besides solving the classical prob- 
lems of modern cosmology (e.g. the flatness, entropy, 
and horizon problems), it provides primordial perturba- 
tions in perfect agreement with the observations @, 9 ■ 
A challenge for inflationary models is to protect the re- 
quired flatness of the inflaton potential against radiative 
corrections. This protection can be provided by an ap- 
proximate shift symmetry as in models of natural infla- 
tion (J,[5|. The symmetry can be broken by a controllably 
small amount, the most known example of this is the case 
of an axion field acquiring a potential from instantons. 
The application of this to inflation, however, requires a 
greater than Planckian axion decay constant / [(| , which 
may not be stable against gravitational corrections Q 
and which may be impossible to realize in string theory 
8[ . Proposed solutions to this problem include using two 
j| or more [lOj axions, which provide an effective large 
scale evolution even if the decay constants of the original 
axions are sub-Planckian, requiring nontrivial compacti- 
fications in string theory [ll| , suitably coupling the axion 
to a 4- form modifying the axion kinetic term [l3j], 
and slowing down the axion evolution through particle 
production [14l |15| as in warm inflation [16| . 

In particular, in the mechanism of [14] the dissipation 
occurs through the production of a U(l) field coupled to 
the inflaton \ through the interaction x^F (where F is 



the U(l) field strength, and F its dual). Ref. [IjJ showed 
that this coupling can also affect the background evolu- 
tion if the U(l) field is replaced by a SU(2) field with 
a nonvanishing vacuum expectation value (vev). Specif- 
ically, due to the interaction with the vev of the vector 
multiplet, the inflaton can be in slow roll even if its poten- 
tial would otherwise (i.e., in absence of this interaction) 
be too steep to give inflation. Such a model has been 
dubbed chromo-natural inflation in [17]. Q 
The model is characterized by the action 



S = 



d 4 x. 
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where x is the axion inflaton and F® = d^A" — d v A a ^ 
ga ahc A\A 



We use the convention e 0123 = 1 for the 
Levi-Civita tensor. The vector field has the vev0 



A a 
^0 



0, A1 = 8?a(t)Q(t). 



(2) 



This mechanism has been extended to a Chern-Simons interac- 
tion in |18|. For recent reviews of vector fields in inflation, see 
0. 

Chromo-natural inflation has trajectories in common with the so 
called gauge-flation model l20t |21| - a model characterized by a 
SU(2) field with a (F^F' 11 ') 2 term besides the usual F^F^ 
Yang-Mills term - in the limit in which the axion is close to the 
bottom of its potential [22, 23]. Perturbations of gauge-flation 
were studied in the last of 1191] ■ 

We follow the standard convention of using greek letters for 
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which is chosen to give isotropic expansion. We note that, 
for a generic theory with vector fields, an isotropic solu- 
tion may be unstable against anisotropic perturbations. 
For instance, this can be the case when a dilaton-like cou- 
pling / (ip) F 2 of a scalar inflaton ip to a SU(2) field with 
the vev ([2]) or to an orthogonal U(l) triplet is arranged so 
to produce a scale invariant spectrum for the vectors. 
Isotropy is instead preserved if the vector fields are mas- 
sive [27| , as in the current model (if the mass arises from 
an explicit breaking, one should also check that the the- 
ory has no ghosts |28|). 

Ref. [ItJ performs a thorough analysis of the back- 
ground evolution of the model. The study of the chromo- 
natural inflation theory at the perturbative level was re- 
cently performed in [291 ] . for 



2 _ 



2g 2 Q 2 > H 2 



(3) 



where m g is the mass of the vector field fluctuations 
in this limit [29[ (as we also discuss below). When 
this condition is realized, the vector field can be inte- 
grated out while, at the same time, leaving its mark 
on the inflationary dynamics (this is an explicit realiza- 
tion of the so called gelaton [3(3] mechanism). In this 
limit, chromo-natural inflation is equivalent to a single 
scalar field P(X, x) = X + X 2 /A 4 - V(x) theory (where 
X = - {dxf /2 and A = 8/V/A 4 ). The non-canonical 
kinetic term precisely encodes the effect of the gauge 
fields. 

In |29j it is also shown that an effective field the- 
ory equivalence of chromo-natural inflation to a non- 
canonical P(X, x) Lagrangian holds as long as m 2 > 
8H 2 . Beyond this limit, a general perturbative study 
of the dynamics of the gauge field and the axion becomes 
necessary in order to test the stability of the theory and 
formulate its predictions. This is precisely the scope of 
the present work. We perform a full linear order study 
of scalar, tensor and vector perturbations and we show 
that the inflationary background solution of the model is 
stable for m g > 2H , and it is otherwise unstable. 

A full phenomenological study of the model is beyond 
our purposes. Nonetheless, we explore the scalar and 
tensor modes production for a given choice of / <C M p 
and A 3> 1 for which the coupling with the vector fields 
is crucial to ensure slow-roll. Specifically, we choose 
/ = 10 _2 M p , and A = 500, while m g /H is a free param- 
eter controlled by the value of g. The stability condition 



m g /H > 2 provides a lower bound for this ratio. Too 
small values are also excluded because they lead to a too 
red spectrum of the scalar modes (we obtained this nu- 
merically; this behavior is also seen in the analytic study 
of [29| in the regime of validity of their analysis). On 
the other hand, the amount of gravity wave signal in- 
creases with m g /H, and a level which can be observed 
in the current or forthcoming experiments 31| can be 
obtained even if the inflaton in the model spans a range 
of some orders of magnitude below M p , evading the so 
called Lyth bound 32] (in contrast to what was expected 
in [13] )■ Therefore, requiring a stable solution, with a 
sufficiently flat scalar power spectrum, and sufficiently 
small tensor modes provides constraints on the parame- 
ters of the model that go in opposite directions. For our 
choice / = 10~ 2 A/p, and A = 500 we could not find any 
acceptable solution. We expect that the situation should 
improve at larger /, where the inflaton potential becomes 
flatter (and the model approaches conventional slow roll 
inflation) . 

The paper is organized as follows: in Sec.|TI]we review 
the background evolution of the model for different (sta- 
ble or unstable) ranges of the theory; in Sec. IHII we start 
the perturbation analysis by introducing the most gen- 
eral decomposition of the metric, the inflaton, and the 
gauge field fluctuations, identifying the physical degrees 
of freedom and verifying that scalar, vector and tensor 
perturbations decouple at linear order; we also quantize 
the system and define the cosmological correlators; in 
Sees. IIV[ fVl and [VTl we study the fluctuations for, respec- 
tively, tensor, vector and scalar modes; finally in Sec. lVIll 
we draw our conclusions. 

We supplement our perturbation analysis in Appendix |A1 
where we show that neglecting scalar metric fluctuations 
does not affect the stability analysis of the model (while 
it considerably simplifies the computations). 



II. THE MODEL AND THE BACKGROUND 
SOLUTION 



Chromo-natural inflation is described by the action (TTJ) 
of a pseudo-scalar field x coupled to an SU(2) gauge 
multiplet. The vector multiplet has the vev The 00 
component of Einstein equations reads 



space-time indices, i,j,k... for space indices, and a, b, c, . . . for 
internal SU(2) indices. The index a should not be confused 
with the scale factor a (t) , that enters in the line element as 
ds 2 = -dt 2 + a 2 {t)S i jdx i dxi = a 2 (r) [-dr 2 + S ij dx i dx^] . We 
denote by dot a derivative with respect to physical time t, and 
by prime a derivative with respect to conformal time r. 
4 This effect originates from the sum of the IR modes that, in 
general, strongly breaks isotropy [24( . It is however possible that 
we live in a realization of inflation where this effect is small, as 
we believe that must be assumed in the computations of |25| . 
Analogous considerations may apply to the model of 



3M 2 H 2 =' i -(Q + H Q) 2 + \g 2 Q± + + V ( X ) . (4) 

Inflation can only occur if the potential V (x) 
is the major contribution to the total energy 
density, i.e. 3H 2 Mp ~ V; the parameters 
e a = Q 2 /M 2 , e 2 = g 2 Q*/H 2 M 2 , e Q = Q 2 /H 2 M 2 
and e x = x 2 /2H 2 Mp are all much smaller than unity. 

The equations of motion for the inflaton and the gauge 



3 



field are: 

X + ZHx - 
Q + 3HQ 



M • X 

7 sm l7 



3^ 2 



Q 2 [Q + HQ\, (5) 



i? 



/ 

2i/ 2N )Q + 2.9 2 Q 3 = ^Q 2 x. (6) 



If we neglect \, Q and H, one can solve Eqs. (0|6|) for x 
and Q 



Q 



g\fQ 2 H - HQ 



fv, x 



3f 2 H 2 + g 2 X 2 Q 4 
HQ (2f 2 H 2 + 2g 2 f 2 Q 2 +g 2 X 2 Q 4 ) 



(7) 



g>~Q 2 fV, x 

3 



3/ 2 # 2 + <? 2 A 2 Q 4 



In these equations we then assume 

3/ 2 iJ 2 < g 2 X 2 Q i , X 2 Q 2 > 2/ 2 . 
and we obtain 



X 



fH {2g 2 Q 3 



gXQ 2 
Q^-HQ 



H 

fV,x 
3gXQ 2 



HQ 



fV X 
gXQ 2 



(8) 
(9) 

(10) 
(11) 



Eq. (fTTj) can be rewritten as an equation of motion for 
the gauge field in terms of an effective potential 



HQ 



dV cS (Q) 
dQ 



H 2 Q 2 



0, V eS 



fHV x 
3gXQ 



where V c s is minimized by 



p 4 sin(x//) 
3gXH 



1/3 



(12) 



(13) 



In the minimum, the two terms of V e g are parametri- 
cally equal to each other. Therefore, the same is true for 
the last two terms in (ITUl) . As a consequence, the first 
term on the right hand side of (JTUJ) dominates over the 
other two in the regime m 2 ^> H 2 studied in [29]. This 
implies that [2!| 



2 g 2 Q 2 > H l 



gXQ 



(14) 



For general values of the parameters, inserting Q n 
Eq. (fTTJ|) becomes 



X 



2 fW ^HW + W'Sj-V^f 3 

3 2 / 3 A 4 / 3 5 2/3 H l/3(_y x )l/3 



(15) 



By looking at Eq. ([5]), one realizes that, when the gauge 
field settles in its minimum, it originates an actual damp- 
ing term for the motion of the axion. Notice also that, 
when Q = Q m im the kinetic energy of the fields can be 
disregarded (see e.g. Eq. (fTTj) ); therefore the slow-roll 
parameter (which would normally receive contributions 



also from cq and e x ) will be, instead, mostly due to ei 
and (.2 



H Q 2 , g 2 Q* 



H 2 ~ M 2 H 2 Mp ' 



Finally we have 

e 2 5 2 Q 4 Q 



V 



He H 2 Mp HMpe 



2Q 



4g 2 Q 3 
H 2 



(16) 



(17) 



From Eq. (|15[) we can compute the number of e-foldings 
■ 4 /3 



N ~ dx 



^i(^) ' 3 2 / 3 A 4 / 3 (sm:r) 1/3 (l + cosx) 2/3 

(mf) A2/3 (1 + cos ^ + 3 2/3 5 4/3 (sin x) 2/3 

(18) 



where we defined x = x/f- 

Let us introduce the parameter y 



y 



A/x 4 



4 \ 2/3 



,3tf 2 M 4 , 

in terms of which the expression (|18[) rewrites 



2 J 2y + z-yz 3 7 \ f 



(19) 



1/3 



(20) 

An upper bound on N can be obtained by "pretending" 
that the slow roll approximation holds at all values of \, 
corresponding to < z < 2 1 / 3 . The resulting upper limit 
is maximized for y ~ 1, where it evaluates to N < 0.6A. 
In the i/« 1 limit, the upper limit acquires the simple 
form N < 

In terms of y, the slow roll solutions read 



H 



V3AL 



vr 



V3Mp 



1/4 / \ 1/3 



^ ■. r V ' I since 



jmfH{x)fA{x) 



A 1 / 2 VVl + cosxy A 1 / 
2/^i 2 y(l + cosa;) 4/3 + (sinj;) 2/3 £T/ 



v^MpA ^(l + cos^'/^sinx) 173 A ^ ^ ' 

(21) 



where we have defined 



1/4 / • \l/3 

(1 + cosx) ' 
2/ (1 + cos a;) 4 / 3 + (sinx) 273 



= VI + cos a; , /a 



•^/y (1 + cosx) 2 ^ 3 (smx) 1 ^ 3 



(22) 



Figure [T] allows to appreciate the accuracy of these slow 
roll approximations. In the figure we show the evolution 
of x f° r f° ur choices of y, comparing the exact evolution 
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FIG. 1. Evolution of x during inflation for four different 
values of y. The other parameters are A = 500 and / = 
10 -2 M p . The evolutions shown correspond to 60 e- folds of 
inflation. 



full numerical background solutions, with analytical so- 
lutions, for which the slow roll approximations are used, 
and we also find excellent agreement. As we discussed 
after eq. (|20p. decreasing y in the y < 1 region, while 
keeping the other parameters (including the initial value 
of x) fixed, decreases the amount of inflation. As a con- 
sequence, y -C 1 can result in sufficient inflation only if 
the inflaton is initially close to the top of the potential. 
Another way to express this is to note that the inflaton 
rolls faster as y decreases (with the other parameters kept 
fixed), as it is clear both from the slow roll expression (|2Tj) 
and from the figure. 

Perhaps the most surprising feature of the model is 
; that, although all the solutions shown in Figure [1] ap- 
— pear to be acceptable inflationary solutions (in all cases 
SOOOOfhe s i QW ro ji so i u tion appears to be an attractor; we note 
however that the background dynamics only probes ho- 
mogeneous departures from the slow roll solution), and 
although they only differ from each other by the value of 
y, the background solution with y = 5 is unstable, while 
the other ones are stable. This emerges from the study 
of the scalar perturbations around these solutions that 
we perform below. 



with the slow roll approximation (|21[), In all cases, we 
take / = O.OIMp, as the main purpose of this model is 
to allow slow roll for / -C M p . This is "compensated" 
by A 3> 1, and in all cases we fixed for definiteness A = 
500. We note that fixing the value of y does specify a 
relation between \i and g, but does not fix these two 
values. As typical in inflationary models, the parameters 
can be specified only from fixing the normalization of the 
power spectrum of the scalar modes to the observed value 
~ 2.5 • 10~ 9 2j. The evolutions shown in the figure 
cover 60 e-folds of inflation. 

We note from the figure that x performs very small os- 
cillations around the slow roll solution. We believe that 
they are due to the fact that Q (t) is tracking a time de- 
pendent minimum (I13|) (we do not show this here, but 
also this tracking is extremely accurate during the slow 
roll phase). The oscillations appear somewhat large in 
the logarithmic scale chosen, but we see that they do 
not lead to any net departure from the slow roll solu- 
tion. Moreover, in Section [VII we compare full numerical 
solutions of the scalar perturbations, obtained using the 



5 The initial conditions for the numerical evolution shown in the 
Figure are chosen as follows: we fix an initial value for \\ we 
then employ (12 1 li to have the initial value for \ an d Q. We set 
Q = 0, and we then obtain the initial value for d by solving the 00 
Einstein equations exactly at the initial time. We make sure that 
the initial value of \ leads to more than 60 e-folds of inflation 
(so that the slow-roll solution can be achieved; we note that the 
system starts slightly displaced from the slow roll solution, since 
we set Q = 0; the displacement is however very small, since 
Q <C HQ in the slow roll solution, and the background evolution 
quickly approaches the slow roll solution) . The evolutions shown 
in Figure [T] are restricted to the final 60 e-folds of inflation. 



III. LINEARIZED EQUATIONS FOR THE 
PERTURBATIONS 



In this Section we discuss at the formal level how we 
compute, we quantize, and we solve the linearized the- 
ory for the perturbations of the model around the back- 
ground solution discussed in the previous Section. The 
discussion is divided in two Subsections. In the first Sub- 
section we discuss how the perturbations can be divided 
into three groups, decoupled from each other at the lin- 
earized level. In the second Subsection we give the form 
of the quadratic action and we discuss how we compute 
the corresponding linearized theory for the perturbations. 



A. Decomposition 

There are 23 perturbations in the system, one of the 
inflaton, 12 of the SU(2) vector field, and 10 of the metric, 
that we decompose as: 



X 

AH 

A? = 



a (Y a + d a Y) 

(Q + SQ) 8 ai + d t (M„ + d a M) 

+e iac (U c + d c U) + U 



.goo = -a 2 (1 - 20) 
g 0i = a 2 (Bi + 8iB) 

ga = a 2 [(1 + 2V>) Sij + WtdjE + d i E j + djEi + h K 



(23) 
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In this expression, a = 1,2, 3 is the SU(2) index (we 
also denote by a the scale factor, as there is no ambi- 
guity between it and an SU(2) index), while i = 1,2,3 
ranges over the spatial coordinates. We denote as "tensor 
modes" the perturbations ti a and hij, which we impose 
to be transverse (ft/iy = diU a = d a U a = 0) and traceless 
(ta — hn = 0); due to these properties, the tensor sector 
contains 4 perturbations. We denote as "vector modes" 
the perturbations Y a , M a , U c ,Bi,Ei, which we impose to 
be transverse = • • • = d^Ei = 0); due to this, the 

vector sector contains 10 perturbations. We denote as 
"scalar modes" the remaining 9 perturbations. 

We point out that the terms "tensor/vector/scalar" are 
appropriate for the perturbations of the metric and of 
the inflaton, as they indicate how these modes transform 
under a spatial rotation. We extend this terminology 
also to the perturbations of the vector field, following the 
notation of [2l[ , even if, strictly speaking, these terms are 
inappropriate (given that the SU(2) index has been used 
in the decomposition) . We nonetheless adopt it since the 
fact that the vector vev is diagonal ((Af) cx Sf) plus the 
transversality properties that we have imposed guarantee 
that the tensor / vector / sectors that we have defined 
above remain decoupled from each other at the linearized 
level [U. 

We Fourier transform these modes 



S(t, x) 



d 3 k 



(27T) 



3/2 



(24) 



where S denotes any of the perturbations, and we study 
the theory in Fourier space. In our stability study, we 
solve for the perturbations at the linearized level, and 
therefore we study a mode with a given momentum k 
(modes of different momenta are coupled to each other at 
the nonlinear level). We can actually fix the orientation 
of k along the z— axis without loss of generality. Starting 
from a general direction for k, we rotate the coordinates 
so that k = k k, where k is the unit vector along the 
z— axis, and k > 0. After a general rotation, (Af) is 
no longer proportional to Sf; however we can re-obtain 
(A?) =QSi through a global SU(2) rotation. Therefore, 
we can set k x = k y = without loss of generality. This 
choice simplifies our algebra. 

We need to remove the redundancy associated to gen- 
eral coordinate and SU(2) transformations. Under an 
infinitesimal coordinate transformation with parameter 



(25) 



and we remove the freedom of infinitesimal coordinate 
transformations by setting if> — E — Ei = 0. Consider 
instead an SU(2) transformation with infinitesimal pa- 
rameter a a = e a + d a e (with e a transverse). Under this 
transformation, 



U -> U + gQe , U t -> U t + gQ 



(26) 



and we can fix the SU(2) freedom by setting U = Ui = 0. 
Clearly, also other modes of the metric and of the gauge 



field change under these transformations, and different 
gauge choices can be made. Our choices are motivated by 
the fact that (i) they completely fix the freedom, and (ii) 
they preserve all the Sgo^ and SAq modes. These pertur- 
bations are nondynamical, as they enter in the quadratic 
action of the perturbations without time derivative, and 
can be immediately integrated out. 

With our gauge choices, and with k — k z , the decom- 
position (|23p acquires the explicit form 



X = X + Sx 
A\ = a <y u Q + SQ + <+, +t x , d z Mi) 

A 2 =a (Y 2 , t x ,Q + 6Q-t+, d z M 2 ) 



A 3 =a (d z Y, 0, 0, Q + 5Q + d z d z M) 



and 



9»v 



-1 + 20 B x 
l + h. 



We verified explicitly that 
(S X , Y, SQ, M, 0, B), the 



B 2 d z B 
h x 
l-h + 
1 

the scalar 
vector 



(27) 



(28) 



modes 
modes 



(11,2, Mi t 2, -61,2) and the tensor modes t > 



h+, h 



are decoupled from each other at the linearized level. 
Namely, the quadratic action for the perturbations splits 
into three decoupled parts 



S-2 — S\ 



2, scalar 



2, vector 



S: 



2, tensor ■ 



(29) 



B. Quantization of coupled systems and correlators 

We have seen that the total action for the perturba- 
tions splits in a sum of three decoupled quadratic actions. 
Let us denote by y the vector formed by the perturba- 
tions in one of these three systems. We can perform a 
transformation 



J'- = Mij Aj 
so that the action for the array A is of the type 



(30) 



drd A k 



A't A' + A'^A - A^A' - A^ 2 A 



Hermitianity of the action implies that K is an anti- 
Hermitian matrix, and Q 2 a Hermitian matrix; the ma- 
trices obtained in the current model are actually real. 
In the following we quantize the system (f3~Tj). so to ob- 
tain the initial conditions for the modes, and we give an 
expression for the correlators between the modes. This 
discussion summarizes the one of 1331 . We first "rotate" 



V> = RA, 



(32) 



where R is a unitary matrix (so that the kinetic term in 
(|3ip remains canonical), satisfying 



R' = RK. 



(33) 
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As K is real, R can also be taken real, and, therefore, 
orthogonal. We note that R is not uniquely determined 
by this condition, and we fix it by setting R = 1 at the 
initial time T; n ; the goal of the present discussion is also 
to understand when this initial time can be set. As we 
will see, the explicit solution for R is not needed. In 
terms of the vector ip, the action becomes 



1 



drd A k 



fl 2 = R(fl 2 + K T K) R T . (34) 
We then introduce the orthogonal matrix C satisfying 



C T n 2 C = diag (uf, wf, u\ 



and we decompose 



(35) 



(36) 



where Ot/oj destroys / creates a quantum with the fre- 
quency u>i. These operators satisfy the algebra 

a i (k S j,a](p)]=6^(k-f)s ij . (37) 

For the systems that we study, in the deep sub-horizon 
regime (aH) in <C k 



Qf n = (n 2 + K T K). = fc 2 1 + a 2 H 2 C 



(38) 



with C constant at leading order in slow roll. Therefore, 
provided that (lf n ~ A: 2 1, we can set the initial conditions 
in the adiabatic vacuum 



-i f T drui ~i II drui 

h~ ^^U = - 



2uj 



2u 



(39) 



where U is a constant arbitrary orthogonal matrix (at 
the practical level, Eq. (|39|l is an approximate solution 
of the equations of motion at early times; the closer f2 2 n 
is to fc 2 l, the better this approximation is. This deter- 
mines how early Ti n needs to be taken). The matrix U 
is unphysical, as it drops from the equations of motion 
for the modes, and from the physical correlators (which, 
as we will see, are given in terms of hh^). The freedom 
associated to U is the generalization to N fields of the 
freedom of changing by a constant phase the wave func- 
tion in the single field case. In the final expression in 
(f3"9"| we have used the freedom associated to U to set the 
wave functions to be real at the initial time. 
Combining (|32|) and (f36|) . we have 



V = R T Ch. 



(40) 



Using the fact that Clf n 
C ~ 1, we arrive at 



fc 2 l at T; n , so that initially 



1 



2k 



(41) 



We start from these initial condition and evolve the equa- 
tions of motion for the modes following from pip 



v" + 2KV + (n 2 + K')V = 0. 



(42) 



In this discussion, all the expressions are given in 
Fourier space. Let us denote by Yi the original fields 
in real space, 



Yi 



d 3 k 



_ik-x 



(27T) 



3/2 



(43) 



Wc have the two point correlation functions 

Cij (x, y) = \(Yi (t, x) Yi (r, y) + Yj (r, y) Y t (r, y) ) 



r dk 
J T 



dk sin (kr) 



k 1,t ' ' J 
where the power spectra are given by 



r=\x-y\ 



(44) 



2tt 2 
2tt 2 



Re 



(MW*M" J 



(45) 



These correlators are the theoretical prediction, to be 
confronted with the statistical average of the correspond- 
ing quantities. The power spectra (j45"T) are the general- 
ization to a system of N fields Yi of the standard power 
spectrum of single field inflation. 

To summarize, starting from the original fields yi in 
momentum space, we perform (|30|) to have a canonical 
kinetic term for Aj. We then decompose this field in 
terms of annihilation / creation operators of the phys- 
ical quanta in the system (the particles of frequencies 
uji). We work in terms of the coefficients of this de- 
composition. We can set the initial conditions (|4Tj) for 
these coefficients, provided that the initial time is chosen 
sufficiently early such that fi 2 + K T K ~ k 2 . Starting 
from these initial conditions, the coefficients T>ij evolve 
according to (|42|) . We will see that for some choice of pa- 
rameters some of the perturbations become exponentially 
large on a timescale <C H^ 1 , signaling an instability of 
the linearized theory. For values of parameters leading 
to stable solutions, the coefficients enter in observable 
quantities through the power spectra (|45[) . 



IV. TENSOR MODES 



We introduce the two doublets 



■ iA, 



A 



iA x 



~ u ~ V2 ' " R " V2 
for the left and right helicities, where 



V2at + 



A, 




(46) 



(47) 
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The action for the tensor modes splits into two separate 
actions, one for the left and one for the right helicity 
doublet, which are formally of the type ([31]) . The action 
for the left-helicity doublet is characterized by 



HQ) 



" //(), 7x 



(48) 



and 



n{ 1 = k 2 -2U 2 - —.{Q' + UQ) 



2 , 3g 2 a 2 Q 4 



^1 1 2 — (xk 



k 2 + a 2 H 2 
2gQ 2 



M 2 



-2 + 



M 2 



'2 



2M 2 



JM, 



" H0 1 7x 



1 



Q 22 = k 2 — aki 2gQ + — ax' ) + g'-jQo-x' 



A 

T 



A 

•r 



k' z - akHA + a 2 H 2 B 



(49) 



where the first expression for each coefficient is exact, 
while the second one is obtained from the slow roll ap- 
proximation (|2ip (we note that Q <C HQ). The quan- 
tities A and B are both of O (A ) . Using the slow roll 
result (fTT|) with Q given by (fT3]). they evaluate to 



A: 



AgQ 
H 



H 2 



2g 2 Q 2 



2g 2 Q 2 
H 2 



B 



H 2 
9 2 Q 2 



(50) 

The action for the right-helicity doublet is related to (|4T)]) 
by fc — ¥ —k, signaling the breaking of parity invariance 
induced by the evolution of the pseudo-scalar inflaton. 

The interactions between the gravity wave and vector 
field tensor perturbations are slow-roll suppressed. How- 
ever, the effective frequency squared f2| 2 01 *z turns neg- 
ative for an intermediate interval of time next to horizon 
crossing. This leads to a tachyonic growth of ti, in 
this interval of time, and, correspondingly, to a growth 
of /il- The same growth does not occur in the right- 
helicity sector due to the opposite sign of the linear term 
in k. The situation is analogous to that first studied in 
[lij . where the interaction %-F-F between a vector field 
and a pseudo-scalar rolling inflaton results in a tachyonic 
growth of the vector modes of a given helicity. 

We note from ((511)) that the length of the tachyonic 
region increases with increasing 



V2gQ 



H 



H 



sin 1/3 x 



V (1 + cos a;) 



2/3 



(51) 



(the last expression is the slow roll approximation) . This 
corresponds to a larger tensor mode production with 
growing ^f. This corresponds to decreasing y in the 
numerical examples that we show below. 



V. VECTOR MODES 

The vector modes Y\ , Y% , B\ , B2 are non-dynamical, 
and can be integrated out. Namely, they enter in 
the quadratic action of the perturbations without time 
derivatives, and therefore their equations of motion are 
algebraic equations in these variables (recall that we are 
in momentum space). When we solve these equations, 
we obtain an expression for the non-dynamical modes in 
terms of the two dynamical modes Mi and M^. We then 
insert this expression back in to the quadratic action for 
the vector modes, and obtain an action for only. In 
other words, the non-dynamical modes do not introduce 
additional degrees of freedom, but are completely deter- 
mined by the dynamical ones. After integrating out the 
non dynamical modes, we define 

Mi =F 1 V 1 + iF 2 V 2 , M 2 = iF x V x + F 2 V 2 (52) 

where 



F 



1,2 



Ml (±k + gaQT + g 2 a 2 Q 2 (M 2 + 2Q 



V2gkM p a 2 Q 



(53) 



(with the upper + sign corresponding to F\ and the lower 
— sign corresponding to F 2 ). The modes Vi t2 are the 
canonical modes of the system, and are decoupled from 
each other at the linearized level: 



1 , , 

S2, vector = - / dnl 'l: 



Wi\ 2 -n 2 \v+ 



WL 



ntJV-\ 



(54) 



An explicit computation gives 



A 



Mp [c 4P 4 ± c 3P 3 + c 2P 2 ± cip + c Q ] 



a 2 P ± f XP+ M 2 p p 2 ± 2gM 2 Qp + 2g 2 Q 2 (M p 2 + Q 2 ) 

where p = — is the physical momentum, and the coeffi- 
cients Ci are functions of background quantities. We note 
that 



n 2 v+ ( P ) = n 2 v _ (-p) 



(56) 



6 At late times, the k— dependence becomes negligible, and the 
mass term for the two vector perturbations t L / R reproduces the 
value m g given in [2S| . 



We also note that the numerator in (j55j) is always pos- 
itive. The exact expressions for Ci are readily obtained, 
but they are not illuminating. We report here the leading 
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expressions in the slow roll A > 1 expansion: 



c 4 = H 2 



C3 



c 2 = H 4 



a = H 5 



2g 2 Q 2 
H 2 

g 2 Q 2 

H 2 

9 2 Q 2 
H 2 

g 4 Q 4 



o 



G.9 



1 

A 

H fH 

, A^ B 2 Q 2 



O 



O t 



i9 H +6 jH 



CO 



A X 



(57) 



Using (HU, it is immediate to see that, in each of the 
above expression, the dominant term in the square paren- 
thesis that we have written explicitly is a O (1) coefficient. 
Moreover, the second term on the right hand side of (f55|) 
evaluates to ±pH O (f). Therefore, at leading order in 
slow roll, 



v± 



p~» H 



Uh -a 2 H 2 



V 



3/4 



p < H 



(58) 



in the deep sub-horizon and super-horizon regime, re- 
spectively. We conclude that the vector sector is stable. 

We note that the super horizon limit of (|58)) is pre- 
cisely the slow roll expression of a 2 m 2 g in the m g ^> H 
limit. The two dynamical vector modes originate from 
the perturbations of the vector field, and the expression 
for the mass in this limit coincides with what found in 
0. 



VI. SCALAR MODES 

After gauge fixing, the system of scalar perturbations 
comprises of 1 mode from the inflaton, 8\, 3 from the 
SU(2) field, Y,SQ,M, and 2 from the metric, <j) and B. 
Among, these 6 modes, Sx, SQ and M are dynamical, 
while the other three modes are nondynamical; namely 
they enter in the quadratic action without time deriva- 
tives, and they can be integrated out (see the discus- 
sion at the start of the previous section). Integrating out 
the metric perturbations makes the algebra extremely in- 
volved; for this reason, in the study presented in the main 
text we make the approximation of setting <p = B = 
from the start. We then integrate out Y, and we de- 
note the resulting quadratic action by Scalar- We also 
performed the full computation, which we present in Ap- 
pendix [A] We denote by S^scaiar-fuii the quadratic ac- 
tion obtained by including all the modes and by inte- 
grating out the 3 nondynamical ones. Both ^2, scalar and 
^.scaiar-fuii are functionals of the 3 dynamical modes. 

We expand all entries in these actions in slow roll. As 
we show in Appendix all entries of the matrices of the 



two actions agree at the leading order in this expansion 
at all scales (namely, for all values of H/p), with a single 
exception. The exception is the the 11 coefficient of the 
frequency matrix Q 2 , for which the agreement is excellent 
only up to H < lOp. The discrepancy that takes place 
afterwards is surely completely irrelevant for the stability 
study that we perform here (as we will see, the instabil- 
ity, if present, manifests itself deeply inside the horizon). 
Not surprisingly, the metric perturbations do not affect 
the stability of the background solution. Very likely, this 
disagreement has also no significant consequence for the 
power spectra that we show below, since it manifests it- 
self only after the power has frozen (see the appendix for 
a detailed study). 

We therefore set <f> = B = 0, and integrate out Y. We 
then define 



<5x 



SM 



Ai 

a 

A 2 

gaQA 2 + y/fc 2 + 2g 2 a 2 Q 2 A 3 
V2gk 2 a 2 Q 



(59) 



In terms of the multiplct A = (Ai,A2,A3) , the 
quadratic action for the scalar modes acquires, up to a 
total derivative, the form (l3~Tj) . We denote as K s and fl 2 
the 3x3 matrices entering in this action. These matrices 
have the following entries 



K s ,u g\Q 2 



a 

K s , 13 
a 

a 



V2f 



g 2 XQ 3 



V2f^ P 2 + 2g 2 Q 2 



= 



(60) 



and 



O 2 



g 2 X 2 p 2 Q 4 



f 2 (p 2 + 2g 2 Q 2 



g 2 Q 4 

M 2 



r 



XX 



2B 1 



Q + HQ 



2M 2 



M 2 



SgXHQ 2 x/2g\QQ 



O 2 

"s,13 



V2f 
x/2X 



i 

g 2 HQ 3 



2vV + 2 5 2 Q 2 



+ 



2p 4 + 3g 2 p 2 Q 2 + 4g 4 Q 4 



2 {p 2 + 2g 2 Q 



2/02 > l 3 / 2 



(q 



HQ, 
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O 2 

"s,22 



o 2 

"s,23 



o 2 



4g 2 Q 2 



gXQx 



Vp 2 + VQ 2 



2.9Q - 



4 5 2 Q 2 (p 2 +ff 2 Q 2 ) 3 Ap 2 Qx 



6g 2 p' 



(q 



2. 9 2 g 2 



/ (p 2 + 2 5 2 Q 2 



+ T~ (61) 

(p 2 + 2g 2 Q 2 ) 2 

We stress that these expressions are obtained by disre- 
garding the scalar metric perturbations, but that they 
are otherwise exact. We can also verify that the eigen- 
values of the {i,k} = {1,2} part of f2|- are ~ a 2 m 2 in 
the super horizon regime and for m g 3> H, in agreement 
with [29[. 

We solve the theory specified by this action following 
the steps outlined in Subsection IIII Bl We assume that, 
after inflation, only the inflaton field provides a sizable 
contribution to reheating (we note that the energy in 
is much smaller than the inflaton energy during infla- 
tion). In this limit, we have the curvature perturbation 

c 



-S-(5x, with the power spectrum 



h 2 v 11 

x 2 



(62) 



where V is given in (14"51) . with being the three fields 
on the left hand side of (|59|) . 

In Figure [21 we then present the time evolution of V\\ 
for a single mode (a given k) and for the same choices 
of background parameters that we used in Figure [1] for 
the background evolution. For definitcness, we consid- 
ered in all cases the mode that leaves the horizon 60 
e-folds before the end of inflation, and we denoted the 
corresponding comoving momentum by k^. We choose 
the initial time of the evolution such that the mode is 
deeply inside the horizon at the start, and the last term 
in ([55)1 is negligible. We observe that the choice y = 5 
leads to an instability of the linearized theory, while the 
other three cases are stable, and are characterized by a 
power that freezes outside the horizon, as in the standard 
inflationary models. 

The instability of the y — 5 choice manifests itself 
at some scale inside the horizon, namely ^unstable = 
k/a uns ta,ue H. Even the stable solutions show a differ- 
ent evolution at this scale (note in the Figure the small 
change in the slope of P (t) for the stable cases at the 
same scale at which the y = 5 solution becomes unsta- 
ble). We can see this analytically, by considering the 
approximated expressions of the scalar system inside the 
horizon. We will obtain the scale AH times a numerical 
factor, where 



A 



A> 1 



(63) 



(in the numerical examples that we have shown, A = 
\/2 10 3 y -1 / 4 ). Let us discuss the approximation in more 
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o 1 e-05 
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FIG. 2. Time evolution of the power (normalized to one at the 
initial time shown), for a mode that leaves the horizon 60 e- 
folds before the end of inflation, and for the same background 
evolutions shown in Figure [T] 



details. First of all, using the slow roll conditions (|21[) . 
we find g 2 Q 2 ~ y- z l 2 f A H 2 < A 2 H 2 in all cases that 
we have studied. Therefore, we can disregard g 2 Q 2 in 
comparison to p 2 inside eqs. (|60| and (|6T|). From this, 
and from the slow roll approximations (|21[) . we obtain, 
in the sub-horizon regime 



K. 



s,12 



gXQ 2 
V2f 



= O (AH) 



H 



= O I AH — } < K 12%s 
\ p 



(64) 



and 



O 2 



O 2 

"s,12 



o 2 



2 X 2 Q 4 



r- 



= O (A 2 # 2 



O (AH 2 

V2\p 
T 



HQ = (AHp) 



O 2 O 2 



O (H 2 ) 



n 2 

"s,23 



= 0(H P ) 



(65) 



We want to solve the evolution equations (|42l) until 
p ~ AH. We can do so by rewriting them as 



-V, 



n 2 ., 



kj 







(66) 



by inserting only the terms written explicitly in (|64j) and 
(|6T)]) . and by treating these terms as constant. All these 
approximations amount in considering only terms that 
contribute to the dynamics at O (p 2 , AHp, A 2 H 2 ). We 
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note for instance that T> = O (pT>) in this regime, so 
that it is consistent to set ft 2 12 = in the approximated 
equation, while retaining if Sj i2 and the dominant term of 
f2 2 13- We also note that the fastest evolving coefficient 
in the matrices is p that evolves on a O (if _1 ) timescale; 
therefore time derivatives of the terms in and (|65|) 
introduce at most terms with an additional factor H, 
which are therefore suppressed in the p > AH and p ~ 
AH regimes. 

Performing these approximations, the system (|66l) re- 
duces to a set of linear second order equations with con- 
stant coefficients. The equations split in three sepa- 
rate groups, one for the complex unknowns T>n (with 
i = 1,2,3), one for the the complex unknowns T>i2, and 
one for the three complex unknowns T>i%. In each group 
we need to solve three second order differential equations, 
and therefore we have 6 possible solutions; we note that 
the three groups have the identical set of equations, and 
they differ only in the initial conditions. Therefore the 
solutions are of the type 



(67) 



where the matrices C a are integration constants (so to 
match the initial conditions) and £ a are constant num- 
bers. The system is unstable if any of the £ a has a real 
and positive part. By solving the system, one can see 
that the only coefficient that can possibly be real is 



3g 2 X 2 Q 4 



2/2 



XQ 



! 2p 2 (H 2 +g 2 Q 2 



9 AVQ6 



1/2 



(68) 



and therefore 



stability <S> p 2 > p 2 = (2H 2 - g 2 Q 2 ) (69) 



where, using the slow roll condition, 



H 2 



= AV /2 /ifiH2- 



3/2 



V 



(70) 



If the expression in parenthesis is negative the back- 
ground solution is stable. Otherwise, there is an insta- 
bility at sufficiently large wavelengths. Therefore 



stability y < 



sin 2/3 x 



2(1 + cos a; 



,4/3 



(71) 



For any fixed y, this condition is violated at sufficiently 
small x. This provides an upper bound to the amount of 
inflation in the model (for any given choice of parame- 
ters). One needs to verify that this upper bound is com- 
parable with the required amount of inflation. For the 



evolutions shown in Figure [TJ the inflaton field, at 60 e- 
folds before the end of inflation, evaluates to \ ~ 2.45/ 
for y = 5, to x - 2 -W for V = 1, to X - 2.00/ for 
y = 0.4, and to x — 1-56/ for y = 0.1. Correspond- 
ingly, the rhs of the condition (fTTj) evaluates to ~ 2.62 
for y = 5, to ~ 1.43 for y = 1, to ~ 0.95 for y = 0.4, 
and to ~ 0.49 for y = 0.1. The criterion (|71[) therefore 
indicates that the choice y = 5 does not lead to a sta- 
ble inflationary solution of 60 c-folds (recall that / and 
A are fixed to 10~ 2 M p and to 500, respectively, in this 
example), while the other choices do. This is in perfect 
agreement with Figure [5J 



Eq. ([69)1 confirms that the instability, if it exists, takes 
place at p < p = AH times an order one factor. For ex- 
ample, for the unstable y = 5 choice, the criterion (|6"9")l 
gives an instability starting at p ~ 2.67A77, correspond- 
ing to aH/p ~ 0.00025, in excellent agreement with the 
evolution seen in the Figure. 

Finally, let us verify that the instability is extremely 
fast. Let us assume that we have a background solution 
with 2H 2 > g 2 Q 2 . At any given moment, modes with 
p > p have Re£ = 0, and are therefore stable. Moreover, 
the instability is also negligible at very large scales, given 
that £ — ¥ for p — ¥ 0. However, for p < p, a quick study 
of eq. (|68p shows that £ = O (AH ) for p smaller than, but 
parametrically equal to p. (for 2H 2 3> g 2 Q 2 , the maxi- 
mum of £ is obtained for p = p/2). This corresponds to 
a very short instability time ~ -Jv <C H~ x . Therefore, 
each mode experiences a strong instability while still in- 
side the horizon, and we conclude that the background 
solution is unstable whenever the condition (1711) is vio- 
lated. We conclude this Section by showing in Figure 
[3] the power spectrum for the stable y = 0.1,0.4,1 
configurations, which we discuss in the next Section. 



le-06 




1e-08 



1000 



k/k„, 



FIG. 3. Power spectra for A = 500 and three different values 
of y. The spectral index (defined as Pq oc k ns ~ 1 ) is n s ~ 0.81 
for y = 1, n s ~ 0.92 for y = 0.4 and n s ~ 0.96 for y = 0.1. 
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VII. CONCLUSIONS 

We performed a complete study of the linear order 
quantum fluctuations for the chromo-natural inflation 
model. We separated the metric, gauge and axion fluc- 
tuations into scalar, vector and tensor modes, verifying 
explicitly that they decouple at the linearized level. We 
computed their equations of motion and worked out the 
quantization of the system. The tensor sector consists 
of the two gravity wave polarizations plus two modes 
from the gauge field. The gauge mode of one helicity 
becomes tachyonic for some finite interval of time next 
to horizon crossing, sourcing one gravity wave helicity. 
The vector sector consists of two dynamical modes that 
originate from the vector field, and that remain perturba- 
tively small. The scalar sector contains three dynamical 
modes, one originating from the inflation and two from 
the gauge field. 

We showed that, for some parameter choice, one of the 
eigenfrequencies of this system, that we denoted by £, can 
become imaginary inside the Hubble horizon, leading to 
a fast instability. Let us compare this with a standard 
result in inflation. Specifically, let us consider the grav- 
ity wave amplitude dgij = a 2 hf } T . The canonical variable 
h c oc ah TT obeys h" c + £l 2 h c = with the dispersion rela- 
tion f2 2 = a 2 (p 2 — 2H 2 ). As a consequence h c oc a out- 
side the horizon, corresponding to a frozen amplitude of 
h TT . Therefore, although the frequency of the canonical 
mode becomes tachyonic outside the horizon, its magni- 
tude is not large enough to lead to a physical instability. 
The situation is analogous for a test scalar field or for the 
scalar perturbations in the standard case. 

For the case at hand, the eigenfrequency £ is given in 
eq. ([55)) . In this discussion, for illustrative purposes, let 
us approximate the full expression obtained in eq. (|68[) 



with £ 2 



V 



cA 2 H 2 , where A > 1 is defined in (l63l) . 



and c an order one factor. No instability appears if c is 
negative. Otherwise, an instability takes place for modes 
of physical momentum p = O (AH) . The timescale of the 
instability is = O (j^jj) <C H~ x , (or, equivalently, 

the amplitude of the scalar modes grows proportionally 
to a large power of the scale factor) which indeed cor- 
responds to a fast instability. This heuristic discussion 
reproduces the results obtained from the precise form of 
£ (see Section[VT]for the precise computation). Moreover, 
we also performed a fully numerical and exact study of 
the scalar perturbations that confirms these analytical 
results. 

The parameters of the model can be chosen so that 
£ remains real inside the horizon (equivalently, c < in 
the heuristic expression) so that the instability is avoided. 
This corresponds to choosing m g = \[2gQ > 2H, where 
Q is the vev of the vector field, and g the SU(2) coupling. 
Quite interestingly, the quantity m g coincides with the 
mass of the fluctuations of the vector field in the m g ^> 
H regime p9j . and it approximates well this mass for 
m g > H. Therefore, the inflationary solution is stable if 
and only if the vector field is sufficiently heavy. 



We have obtained the boundary of the stable region 
m g 3> H analyzed in [29j and our formalism can be read- 
ily employed to study the phenomenology of the bound- 
ary region. This study is beyond the purposes of the 
present work. For illustrative purposes, we have how- 
ever computed the power spectrum of £ for some sample 
choices of parameters. Specifically, Figure [3] shows the 
power spectrum for an axion decay constant / = 0.01Af p , 
for A = 500 (we stress that the main motivation for the 
model is to provide inflation for a sub-Planckian axion 
decay constant, and that this can be obtained for suf- 
ficiently large A), and for three choices of y. We note 
that, among those shown in Figure [3l only the power 
spectrum obtained for y = 0.1 is sufficiently flat to meet 
the observational bounds 0, Q . We numerically found 
that the spectral tilt is a decreasing function of y in the 
0.1 < y < 1 interval (equivalently, it is an increasing 
function of m g /H in this interval). This behavior can 
also be seen in the analytic relation given in [29( in the 
large m g regime. Therefore, both the requirement of sta- 
bility and of flatness of the power spectrum pose a lower 
bound on m g /H. 

On the other hand, the discussion around eq. (|5Tj) 
leads to the conclusion that a large m g /H leads to a de- 
tectable or ruled out gravity wave signal. The enhanced 
signal is parity violating, which results in nonvanishing 
TB and EB correlations. A measure of the net handed- 
ness of the tensor modes is 



|A%I = 



Pl - Pn 



Pl + Pr 



(72) 



where Pl/r is the power spectrum of the left/right- 
helicity gravity wave modes. The corresponding observa- 
tional bounds have been studied in [34l |35H . In Figure 2] 
we compare the bounds presented in [35| with the value 
of r and A% obtained in this model, for our choice of 
/ = O.OlMj, and A = 500, and for different values of y in 
the 0.35 < y < 0.7 range (the tensor power spectra are 
obtained from the 1 1 element of P5")l for the left and right 
helicity tensor sectors, following the procedure outlined 
in Subsection IIII B|) . Greater values of y do not lead to 
a visible gravity wave signal even in a cosmic variance 
limited experiment, while lower values are ruled out by 
the current limit r < 0.13 [111]. 

We point out that all values of y considered in this 
plot are actually ruled out because they lead to a too 
small value of n s (the largest value n s ~ 0.93 is obtained 
for y = 0.35, and n s then further decreases at larger y). 
We stress that our choice of / = 0.01Af p and A = 500 
is only dictated by the requirement of / <C M p , A 3> 1, 
so that the coupling to the gauge field plays a relevant 
role for the inflation dynamics (which, in turns, is the 
initial purpose of the model [13] )■ Obviously, these values 
can be changed, and our phenomenology discussion has 
the only purpose of understanding what kind of limits 
can be imposed on the model. We expect that a viable 
region will exist at larger / (and smaller A), as the model 
becomes closer to a free inflaton model slowly rolling on 
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a flat potential. 
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FIG. 4. Red/solid curve: Value of Ax and r obtained for 
/ = 0.01A/p and A = 500, and for different choices of y in 
the 0.35 < y < 0.7 interval (successive points along the curve 
denote 0.05 increments in y); black/dotted vertical line: r < 
0.13 bound from ||; the other black/dotted curves are the 
la detection lines for the Planck (P), SPIDER (S), CMB-Pol 
(C), and a cosmic- variance limited (CV) experiment. The 
signal needs to be above a line to be detectable at la by that 
experiment. These experimental forecasts are an approximate 
copy of the lines shown in Figure 2 of [35l ]. 



At the theoretical level, it is interesting to note that the 
model can result in an observable gravity wave signal in 
the near future, even if the infiaton spans a range which 
is some order of magnitudes smaller than the Planck 
scale. This evades the Lyth bound (32[, that states that 
r > 0.01 is possible only for an excursion of the infiaton 
of O (M p ) during the last ~ 60 e-folds of inflation. We 
note that the model we have studied evades the bound 
because of the x^F interaction. Quite interestingly, the 
only other examples that we are aware of where r > 0.01 
can be achieved even if the inflation evolution is orders 
of magnitude below the Planck scale Q are those studied 
in refs. [Til . l39j . which are characterized by the same 
pseudo-scalar interaction. As we already mentioned, ref. 
[14| is also characterized by a ~ / <C M p evolution due 
to the damping from gauge field production. In this case, 
an interesting parity- violating gravity wave signal can be 
generated [40j from the produced gauge quanta, but one 
needs to evade the simultaneous generation of too many 
non-gaussian scalar density perturbations [|l| . Ref. (40| 
achieves this by considering a large number of gauge fields 
(this reduces non-gaussianity by the central limit theo- 
rem) or by introducing a curvaton field. Ref. shows 
that r > 0.01 can be obtained, and the non-gaussianity 



7 We note that inflationary potentials leading to r J> 0.01 have 
been constructed where the inflationary range is smaller than, 

but stiii of o (m p ) mm. 



limit can be respected, if the rolling-scalar is not the in- 
fiaton. 

An important difference is that, however, in the mecha- 
nism of [H], Hil the tensor modes are produced at the non- 
linear level by the vector fields produced by the rolling 
infiaton. For the model [13], the production occurs al- 
ready at the linear level, due to its mixing with the vec- 
tor mode induced by the vector vev and the non-abelian 
structure (<?, Q ^ 0). Quite likely, also in this model 
the vector modes can source a significant amount of 
scalar density perturbations at the non-linear level. This 
may reduce r from the level studied here, although a too 
large m g will still likely be ruled out by the significant 
gravity wave production. It is possible that, for regimes 
resulting in acceptable r, the sourced scalar modes will 
lead to interesting levels of non-gaussianity and primor- 
dial black holes as those obtained in pd>l44l |. It is also 
possible that, for some choice of parameters, the vector 
field production will "self-regulate" (with a consequent 
decrease of r) due to its backreaction on the inflationary 
dynamics, as visible in the background evolutions stud- 
ied in [141 l45j ] (a too large production may slow down 
the infiaton, and this may in turns decrease the vector 
production). All these interesting possibilities remain to 
be studied. 

Note added: In the first version of this manuscript, 
we pointed out that (i) the inflationary solution of (l7| 
is stable if and only if m g > 2H, and that (ii) the 
tensor-to-scalar ratio r in this model is enhanced with 
respect to the case of a free infiaton, leading to violation 
of the Lyth's bound [32[ for some choice of parameters. 
Both these claims are confirmed by the present analy- 
sis. The stability study is unchanged with respect to the 
first version. The study of the tensor modes, and the 
phenomenology considerations that follow from it, have 
instead been updated, to include the helicity violating 
terms in the tensor action that were erroneously miss- 
ing in the first version. Such terms result in a further 
increase of r and in a helicity violating gravity wave sig- 
nal. The relevance of these terms was pointed out in [46j . 
that appeared on the archive between the first and the 
current version of this manuscript. Ref. (4(| agrees with 
our limit m g > 2H for the stability of the inflationary 
solution, and our revised analysis agrees with the effect 
of the helicity violating effects found in (46|. 
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Appendix A: Including scalar metric perturbations 



roll suppressed quantities: 



In Section lVTl we studied the system of linear perturba- 
tions for the model, disregarding the perturbations of the 
metric. Here we summarize the results for the full sys- 
tem, and we confirm the accuracy of the approximation 
made in the main text. 

We start from the full set of 6 scalar perturbations (af- 
ter gauge fixing), including the 2 modes from the metric 
perturbations (<j> an d B) that we have (artificially) set 
to zero in the main text. Both these modes are non- 
dynamical, and we integrate them out, together with the 
other non-dynamical perturbation Y . We are left with 
three dynamical modes, and we "rotate" them as in eq. 
(15^1) of the main text. In this way, we obtain the full 
quadratic action for the perturbations. It is formally of 
the type 



<S2, scalar-full = — I C?TC? 3 fc 



A'+Cs.fA 



(Al) 



Namely, the modes are not the exact canonical 
scalar variables. However, as we now show, they provide 
a very good approximation to the canonical modes in the 
slow roll regime. More in general, the matrices of the full 
action (jAll) are extremely involved. We studied them in 
slow roll approximation. Specifically, using the slow roll 
approximation (12T1) , we can cast all the elements of these 
matrices in the form 



(A2) 



where the coefficients Cj and dj only depend on the pa- 
rameters of the model, and on slowly evolving back- 
ground quantities. We computed the leading order ex- 
pression for these coefficients in the slow roll approxima- 
tion. For example, we obtain 



(C s ,f)u 



g 2 Q 4 x 2 /M 4 



2H 2 p 2 + 4g 2 H 2 Q 2 [l + O 



O 



= 1 



( f 2 , ) 

y M 2 y3 /2 X 3 J 



H 1 



(A3) 

This expression is extremely close to one, since, paramet- 
rically, the second term is A~ 3 <C 1 outside the horizon, 
and even more suppressed inside the horizon. In fact, the 
kinetic matrix differs from the identity only up to slow 



(Cs,f) 



12 



O (A 1 / 2 ) H 2 
0(y 3 / 2 X 3 )p 2 + 0{X 3 )H 2 



(C s ,f) 22 

( C 8,f) 23 

(^,1)33 



= 1 



O (A 3 / 2 ) H^JH 2 + O (y 3 / 2 ) p 2 
O (y 3 / 2 A 3 ) p 2 + O (A 3 ) H 2 
O (A) H 2 
0(y 3 / 2 X 3 )p 2 + 0(X 3 )H 2 



O (A 2 ) H^JH 2 + O (y 3 / 2 ) p 2 
O (y 3 / 2 X 3 ) p 2 + (A 3 ) H 2 
O (A 2 ) H 2 
1+ 0(y 3 / 2 A 3 )p2 + 0(A 3 )^ 



(A4) 



Therefore, up to very small slow roll corrections, the 
modes are also the canonical variable of the full 
scalar system. 

Performing the same procedure on K s j , we obtain 



K 



a 



a 



gXQ 2 P 2 + m 2 g [1 + Q(j)] 
V2f P 2 + m 2 g [1 + 0(1)] 



V a 



1 + 



12 
2 2 AQ 3 



l+y 3/2 p 2/ H 2 



[1 + W 



13 



V2f^/p 2 + m 2 P 2 + m 2 [l + 0(I)] 

" / 13 



13 
tf 4 



23 



o 



1 



,1 + 0, 

,9/ 4 A; H 2 \y 15 / 4 X 



1 




where we recall that m 2 = 2g 2 Q 2 . The 12 and 13 entries 
are in excellent agreement with those given in the main 
text. We recall that the 13 entry is much smaller than 
the 12 entry inside the horizon, and it is negligible in the 
stability study. We note that, for the full system, the 23 
element is nonvanishing, while (^ s ) 2 3 = 0. However, this 
element is strongly slow roll suppressed with respect to 
the other two, and completely negligible. 

Proceeding in the same way (for brevity, we omit here 
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the powers of y), we obtain 



M u 1 , EL o(g 2 > 4 - 2 ') 



(^ 2 )l2 



(p 2 + m 2 ) A 



(«2.f) 



(^ 2 ) 



i3 =1 + I ^ 2 E-= o(tf 2 V 



13 



Aj5 2 + m 2 J2Lo°( H2l P 4 ~ 2i ) 



22 



(^s)22 ^ 



1 + 



^23= J , SLoO(^ 2 V 



2^6-2^ 



(^) 



2:) 



(fi 2 ) 33 fc 2 



= 1 



(p 2 + m 2 ) A 

Et o(gy- 2 ') 

(p 2 + ^)aeL o( 5 ¥- 2 ') 

(A6) 



and we see that, for all these entries, the expressions of 
f2 2 f and of fi 2 agree at all scales (namely, for any value 

of H/p) up to subdominant (A -1 ) corrections. For the 
11 entry, we obtain 



s,f 



H 2 (p 2 ) + (H 2 ) , s 

1 + A, J ' ( A7 ) 



(^ 2 ) n — A: 2 p 2 A p 2 + to 2 



Also on this entry, the two matrices are in perfect agree- 
ment during the full sub-horizon regime. However, while 
for the other entries the agreement continues also in 
the super- horizon regime, the 11 entries differ from each 
other for H 2 > O (Ap 2 ) . By evaluating the coefficients 
in (|A7j) . we found that the disagreement starts only at 
H > lOp for all the evolutions studied in the main text. 



To conclude, all the matrix elements, up to one ex- 
ception, of the system of scalar perturbations studied in 
the main text (in which we made the approximation of 
disregarding the scalar metric perturbations) agree at all 
scales (namely, for any value of H/p) with the corre- 
sponding entries of the full system up to subdominant 
terms in a slow-roll expansion. The single exception is 
the 11 entry of fl 2 , for which the agreement persists dur- 
ing the entire sub-horizon regime, at horizon crossing, 
and also up to H < lOp, but not further. This proves 
that the stability study performed in the main text is 
valid also when the metric perturbations are included, 
given that the instability, when present, manifests itself 
deeply inside the horizon. Most likely, this guarantees 
that also the power spectra shown in the main text are 
accurate, since the disagreement manifests itself only af- 
ter the powers have frozen (see Figure [2|) . 
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